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Abstract 

The incoherent pion photoproduction reaction ^d — J- ■K~pp is considered theoretically in a wide 
energy region Eth < E^ < 2700 MeV. The model applied contains the impulse approximation as 
well as the NN- and vrA^-FSI amplitudes. The aim of the paper is to study a reliable way for getting 
the information on elementary 7n— 7>7r~p reaction cross section beyond the impulse approximation 
for 7d— T-vr^pp. For the elementary jN ^ irN , NN^NN, and ttN^ttN amplitudes, the results 
of the GW DAC are used. There are no additional theoretical constraints. The calculated cross 
sections da/dil{'~fd — )• ir^pp) are compared with existing data. The procedure used to extract 
information on the differential cross section da/dQ{'yn^ n^p) on the neutron from the deuteron 
data using the FSI correction factor R is discussed. The calculations for R versus n^p CM angle 
9i of the outgoing pion are performed at different photon-beam energies with kinematical cuts for 
"quasi-free" process 'jn — t- Tr~p. The results show a sizeable FSI effect R ^ 1 from S-wave part 
of pp-FSl at small angles close to 9i 0: this region narrows as the photon energy increases. At 
larger angles, the effect is small — 1| <^ 1) and agrees with estimations of FSI in the Glauber 
approach. 

PACS numbers: 13.60.Le, 21.45.Bc, 24.10.Eq, 25.20.Lj 
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I. INTRODUCTION 

The A^* family of nucleon resonances has many well-estabhshed members several of 
which exhibit overlapping resonances with very similar masses and widths, but with dif- 
ferent spin-parity values. Apart from the iV(1535)l/2^ state, the known proton and 
neutron photo-decay amplitudes have been determined from analyses of single-pion pho- 
toproduction. The present work studies the region from threshold to the upper limit of 
the SAID analyses, which is W = 2.5 GeV. There are two closely spaced states above the 
A(1232)3/2+: A^(1520)3/2- and iV(1535)l/2-. Up to a CM energy oiW ^ 1800 MeV, this 
region also encompasses a sequence of six overlapping states: A^(1650)l/2^, A^(1675)5/2~, 
Ar(1680)5/2+, Ar(1700)3/2-, Ar(1710)l/2+, and A^(1720)3/2+. 

One critical issue in the study of meson photoproduction on the nucleon comes from 
isospin. While isospin can change at the photon vertex, it must be conserved at the final 
hadronic vertex. Only with good data on both proton and neutron targets can one hope to 



disentangle the isosca^ 
resonances (see, Refs. 



ar and isovector electromagnetic couplings of the various N* and A* 
3j), as well as the isospin properties of the non-resonant background 
amplitudes. The lack of 7n — t- vr^p and 7r°n data does not allow us to be as confident about 
the determination of neutron couplings relative to those of the proton. Some of the A^* 
baryons (A^(1675)5/2^, for instance) have stronger electromagnetic couplings to the neutron 
relative to the proton, but the parameters are very uncertain [l| . Data on the — )■ vrA^ 
reactions are needed to improve the amplitudes and expand them to higher energies. 

Incoherent pion photoproduction on the deuteron is interesting in various aspects of 
nuclear physics, and particularly provides information on the elementary reaction on the 
neutron, i.e., '-fn-^irN. Final-state-interaction (FSI) plays a critical role in the state-of-the- 
art analysis of the 'jN — )• vr interaction as extracted from 'yd — )■ ttNN data. The FSI was first 
considered in Refs. {4, 5| as responsible for the near threshold enhancement (Migdal- Watson 
effect) in the A^A^-mass spectrum of the meson production reaction A^A^ — NNx. In Ref. Q], 
the FSI amplitude was studied in detail. Calculations of A^A^- and vrA^-FSI for the reactions 
'jd — !■ irNN can be traced back to Refs. 7-^. In Refs. 9], the elementary 7A^ — )■ vrA^ 
amplitude, constructed in Ref. {tI from the Born terms and A(1232)3/2+ contribution, was 
used in 'jd-^TcNN calculations with FSI terms taken into account. Good descriptions of the 
available deuteron data for charged pion photoproduction in the threshold and A(1232)3/2"'" 
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regions were obtained. 

Further developments of this topic (see 10Nl7l| and references therein) included improve- 
ments of the elementary — )■ ixN amplitude, predictions for the unpolarized and polarized 



(polarized beam, target or both, see 10|, ll2Nl6| and references therein) observables in the 
'yd^nNN reactions, and comparison with new data. Different models for 'yN ^ttN ampli- 



tude were used in the above mentioned papers, i.e., MAID IS] (Refs. 12yi3|), SAID 
(Refs. 13, 15 1), and MAID Q (Ref. 15|). As discussed in Refs. 13 



19| 



ISH. the main un- 



certainties of 'jd — )■ ttNN calculations stem from the model dependence of the — )■ vrA^ 
amplitude. In the latest SAID lol and MAID [2^ analyses, the models for — )■ ttN am- 
phtudes are developed for the photon energies < 2.7 GeV [l9| and < 1.65 GeV |20| . 
respectively. Summary results from the existing jd — > ttNN calculations show that FSI 
effects significantly reduce the differential cross section for n^pn channel, mainly due to the 
pn rescattering, and contribute much less in the charged- pion case, i.e., in n^nn and 7r~pp 
channels. 

The role of FSI depends on the kinematical region considered. In Ref. 2l|, a narrow 
enhancement in the pp-mass spectrum observed in the reaction pp^ppir' with backward 
outgoing 7T^ was explained by the pp-FSl. The result was shown to be model-independent, 
determined only by pp-scattering parameters for the pp pair produced at high momentum 



transfer. In the same approach, it was shown [22| that the observed energy behavior of 
the total cross section of the reaction pp — )■ pprj in the near threshold region can be also 
explained by pp-FSI. In Ref. the meson photoproducton on deuteron was considered 
at high energies {E^ ~ several GeV) and high momentum transferred to final meson. This 
work was focused mainly on special kinematical regions close to the logarithmic singular- 
ities of the triangle NN- and vrN-FSI amplitudes, the latter are strongly enhanced. These 
configurations where the FSI amplitudes dominates may be interesting, say, in connection 
with color transparency hypothesis 23|]. On the other hand, to extract the neutron data, 
we are interested in the opposite case, i.e., when FSI is suppressed. 



In this paper, the role of 



analysis addresses the data 



SI in the 'yd ^ it pp reaction is under consideration. Our 



25| that come from the 76?— t-tt pp experiment at JLab using 



CLAS for a wide range of photon-beam energies up to about 3.5 GeV. The calculated FSI 
corrections for this reaction are further used to extract the 7n — > n^p data that constrain 
the 7A^ — )■ vrA^ amplitude used in PWA and coupled channel technologies. 
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FIG. 1: Feynman diagrams for the leading components of the "fd — )• 'K~pp amphtude. (a) Impulse 
approximation, (b) pp-FSI, and (c) vrN-FSI. Filled black circles show FSI vertices. Wavy, dashed, 
solid, and double lines correspond to the photons, pions, nucleons, and deuterons, respectively. 

In our approach, the 7(i — )■ TT~pp amplitude has three leading terms, represented by the 
diagrams in Fig. [T] impulse approximation (lA) [Fig. [T]^a)], pp-FSI [Fig. [T](b)], and vrN-FSI 
[Fig. [11(c)] contributions. lA and vrN diagrams [Figs. [It^a),(c)] include also the cross-terms 
between outgoing protons. It is convenient to study the FSI effects in terms of the ratio 

Rfsi = {da/dn^p)/{da'^/dn^p), (1) 

i.e., the ratio of the differential cross sections da/dQ^j^p including the full calculations of dia- 
grams [Figs. [Tl^a)-(c)] to the {da^^/dQT,p), associated with lA diagram [Fig. [U^a)], where fi^p 
is the solid angle of the relative motion in the final np system. The ratio Rpsi © depends 
on different kinematical variables. It can be used to extract the differential cross sections 
da/dQ for the reaction 7n — > n^p from the 'yd — > n^pp data. We use the recent GW pion 
photoproduction multipoles to constrain the amplitude for the impulse approximation {26 1 
with no additional theoretical input. While for the pp-FSI and vrN-FSI, we include the GW 



NN [27] and GW vrN amplitudes [28|, respectively, for the deuteron description, we use the 



wave function of the CD-Bonn potential 29|] with S- and D-wave components included. 

This paper is organized as follows. In Section [Tll we describe the model. In Subsec- 
tions III Al and III Bl we introduce the notations and write out the impulse approximation 
terms of the '-fd^TcNN amplitude. In Subsections III CI and III Dl we derive the NN-FSl and 
ttA^-FSI terms of the reaction amplitude, respectively. 

The results are presented in Section IIIII In Subsection IIII Al we compare our numerical 
results for the cross section da/dQ{'yd^TT^pp) with the DESY data and discuss the contri- 
butions from different amplitudes. In Subsection IIII Bt we discuss the procedure to extract 



the cross section da/dQ{'yn-^'n'~p) for the neutron from the 'yd^7i~pp data and define the 
correction factor R. In Subsection IIII C\ we present the numerical results for the factor R 
and discuss the role of the S'-wave pp-FSl. In Subsection IIIID[ we estimate the R factor in 
the Glauber approach. The conclusion is given in Section IIVI 



II. MODEL FOR ■yd Tv'pp AMPLITUDE 



A. Kinematical notations 



Hereafter m, /i, and are the proton, pion, and deuteron masses, respectively; q = 
{E-y,q), pd = {Eci,Pd), k = (w, fc), and pi = {Ei,Pj) {i = 1,2) are the 4-momenta of the 
initial photon, deuteron and final pion, nucleons, respectively; k' = {u', k'), p' = {E',p'), and 
p'- = {E'', p9 the 4-momenta of the intermediate particles. The 4-momenta are shown 
in Fig. [H The total energies E^, Ed, ■ ■ ■ E[ and 3-momenta q, p^, ■ ■ ■ p ■ are given in the 
laboratory system (LS), i.e., in the deuteron rest frame, where = and Ed = md. 

The cross-section element daljd^ii'pp), according to the usual conventions for invariant 
amplitudes and phase spaces (see Appendix IV Ap . can be written in the form 



1 \M^d? , , d'p2 ^ , k^dQi 

""^ = o "^3, dT3 = dT2, dT2 = . (2) 

2 AE^nid {2'nY2E2 [A-n)^ Wi 



Here: M^^ is the yd^n~pp invariant amplitude; jM^^P is the square |M^(i| , calculated for 
unpolarized particles; dr^ is the ttNN phase space element, written in terms of the vrpi-pair 
phase space element dT2 and 3-momentum P2 of the 2nd proton; the factor ^ in da ([2]) takes 
into account that the final protons are identical; ki and Qi are the relative momentum and 
solid angle of relative motion in the irpi system, respectively; Wi is the effective mass of the 
vrpi system. 



B. Impulse-approximation amplitudes 



Let us use the formalism of Ref. [30|], which is similar to that of Gross [31] in the 
case of small nucleon momenta |pp/m ^ m in the deuteron vertex. Then, the impulse- 
approximation term Ma [Fig. (U^a)] of the yd^irNN amplitude can be written in the form 

M„ = M«+Mf, (3) 



Here: Ui is the bispinor (isospinor also) of the z-th final nucleon, uu = 2m; u'^ = T2UcU^ = 

fl 2) — ^ fl 2) 

^212^*, where Uc = 7270 is the charge-conjugation matrix; = ■Ui^2^-y7v " 

plitude of subprocess 'yN ^TrNi^2, and u is the bispinor (isospinor also) of the intermediate 
nucleon with 4-momentum p' = Pd —P2,i] Gn{p') = {^'+ m)/(p'^— m^+iO) is the nucleon 
propagator, where = p^7^; Td{p2,i —p') is the dNN-vertex related to the deuteron wave 
function (DWF) as given in Appendix IV B[ The amplitude Ma is antisymmetric with respect 
to the nucleon permutations in accordance with the Pauli principle. 

Further, we retain only the positive-energy part of the nucleon propagator Gn{p') and 
apply the connection between f ^ and DWF Then, for a given spin and isospin states of 
the particles, we obtain 

= 2yM 5^(7r,mi,ri|M^;^| A,m',r')(m',r',m2,r2|^rf(P2)l"^d), (4) 

m' ,t' 

and the 2nd term is Ma^^ = —Ma^^ (with permutation of the variables of the final nucleons). 
Here: mi, 2, rn', A, and m^ are spin states of the final nucleons, virtual nucleon, photon, and 
deuteron, respectively; vr, ri,2, and r' are isospin states of pion, final nucleons, and virtual 
nucleon, respectively. Substituting isospin states for the reaction jd^n'pp, one gets 

= 2^ J2 [ ("^1 1 I K m') (m', m2 1 ^^(Ps) I ^d) (5) 

- (m2 1 M^^) I A, m') (m', mi | ^aiPi) I "^d) ] , 

where now are the 777, — )■ vr pi amplitudes. The expressions for 

DWF {nil, 1712 \ '^d{p)\^d) are given in Appendix IV B[ The 7A^ — vrA^ amplitudes M^at 
can be expressed through the Chew-Goldberger-Low-Nambu (CGLN) amplitudes [32[ (see 
Appendix IV CI) . The CGLN amplitudes as functions of the npi invariant masses Wi depend 
on the virtual nucleon momentum p' through the relation = {q + p'Y = + Pi)^- Thus, 
the Fermi- motion is taken into account in the 'jd ^ n^pp amplitude ([5]). The matrix 
elements (mi| My at] A,m') are given in Appendix IV Dl 
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C. NN final state interaction 



The NN-FSI term Mf, [Fig. [U^b)] of the 'jd^irNN amphtude can be written in the form 



Mb = —I I - — -- > ■ . 6 

Here: m[ and m'g are spin states of the intermediate nucleons; the notations for mi, 1712, A, 
and 171^ are the same as in Eqs. (for short, we omit isospin indices); is the amphtude 
of subprocess jd^TrNN in impulse approximation 

(m;,m^|M^^| A,mrf) = 2v^ ^ (m;| M^^| A, m') (m', m'al m^), (7) 

m' 

where Mat at is the A^A^-scattering amphtude. The integral over the energy in Eq. ([6]) can 
be related to the residue at the nucleon (momentum V^le with positive energy. Let us 
rewrite the 3-dimensional integral J dpl^ in the AW center-of-mass system. Then, we get 
p[^-m^+iQ = 2W{E- E'+ iO), where W is the AW-system effective mass, E = W/2 = 
^/p^Tm^, pj^ = \pj^\, E' = ^/p^'+m?, p'j^ = \p'j^\, and PAr(Pjv) is the relative 3-momentum 
in the final (intermediate) NN state. We obtain 

M.= /-M^ 1:! . (8) 

J (2vr)3 4:E' W {E' - E - iO) 
(■•■)= ^ {m[,m2\M!^j\X,md){mi,m2\MNN\m[,m2). 

One can rewrite as 

Mh = Mr + M°^ f -- 



dp'N {■■■) 



inS{E'~E) + P—^ 



E'-E 



(9) 



(27r)3 AE'W 

Here, M^" and M^^^ are the contributions from the 1st and 2nd terms, respectively, in 
square brackets in the r.h.s. of Eq. (|9]), where P means the principal part of the integral. 
The amplitudes M^" and M'^^^ correspond to the on-shell and off-shell intermediate nucleons, 
respectively. For M^", we get 

M^^ = i'K 5iE'-E) = -^^- fdn' (■■■), (10) 

where dfi' = dz'dip' {z' = cos 6') is the element of solid angle of relative motion of the 
intermediate nucleons. Consider the 2nd term M^-^^ . Let us use Eqs. (|35|) of Appendix IV Bl 



for ^(^(^2) in Eq. ([7]) and represent the integrand (■ ■ ■ ) in Eq. (jS]) as the sum of two terms, 
proportional to S- and D-wave components of DWF, i.e., ^(^2) and wij)'^). Then, we obtain 

{■■■)= Au{p'^) + Bw{p',) (p'2 = |p'2|), 
A = 2^/m {m\\M^N\\,m'){m\m'2\Su\rnd),{^ ^ ^ 



and B is given by the expression for A after the replacement Su^Sw, where Su and 5*^ are 
given in Eqs. (!35|) of Appendix IV B[ The factors A and B contain and AW amplitudes, 
spin structure of DWF, and depend on the momenta of the particles in Fig. [Il^b). Note that 
the A7V-FSI amplitude ([S]) takes into account the Fermi-motion, since the amplitudes 
My AT and Mat at depend on the intermediate momenta p' and pg, respectively. In the integral 
Jdp'j^ = j dO! dp'j^p'j!} ([9]), we take out of subintegral ^dp'j^ the factors A and B ffTTj) at 
p'n — Pn^ calculate A and B as well as the amplitudes M^n and Matat with the 

on-shell intermediate nucleons. This approximation means that we neglect the off-shell 
dependence of the 7A^ and A7V amplitudes in comparison with sharp momentum dependence 
of DWF. Then, we get 

M^/f = tA ^ =^__ [dQ'(AI^ + BlJ, (12) 

T _ A. (^p'nPn u{p2) J _ J^dp'j^p^ w{p'2. 



nE' E'-E' J nE' E'-E' 

where <f denotes the principal part of the integral. We also include the formf actor f{p'j^) [l 

to parametrize the off-shell ^Sq partial amplitude of pp-scattering and define the integrals 

r(0) _ / dPNPNU{p'2)f{p'N) , , . _ Pn + f^^ 

~f^' W^^' f^P^)-^^:^ ^^^^ 

with /3 = 1.2 fm"^ Q; I^^ = I^°\u{p'2) ^wip'^)). Let us write the terms A and B as 

A = Ao + A^, B = Bo + B,, (14) 

where Aq (Ai) is given by Eq. ( ITT]) when only ^Sq part is saved (excluded) in the pp-scattering 
amplitude Mnn (for Bqi the substitution Su in Eq. ffTTl) is implied). Combining 

Eqs. ([9D-([I3D, we obtain 

The integrals J^, Iu \ Iw \ and / d^l' (fT5|) are carried out numerically. The ATV-scattering 
amplitude is described in Appendix IV El 
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D. vrA^ final state interaction 



The ttA^-FSI term Mc [Fig. (U^c)] of the 'jd^nNN amphtude can be written in the form 
= M« + Mf, Mi^) = - , (16) 

(■■■)= ^ (7^^rl,ml,r2,m2| M^;t| A,md) (7r,r2,m2| M^^l 7r',r2,m2), 

Tr',T^,m'2 

where the integral over the energy is also related to the residue at the nucleon pole (mo- 
mentum P2) as in Eq. ([7j). Here: and are spin and isospin states of the intermediate 
nucleon with 4- momentum is isospin states of intermediate pion; the notations 

ri2, vr, A, and are given above (see Eq. dlD); = (vr, r2, m'2\ M^^| vr', Xg, m'g) is the 

ttA^ — )■ 7rAi'2 amplitude; ^2 is the relative 3-momentum in the intermediate vrA^ system. The 
2nd term Mf^ = —Mc^^ (with permutation of the final nucleons). Substituting isospin 
states for the reaction 7(i— >7r"pp, and making use of Eq. ([7]), we get the integrand (■ ■ ■ ) in 
Eq. fll6p in the form 

{■■■) = 2yM [{mi\M^^\in^7r-p)\X,m') {m2\Ml^}p\m'2) (17) 

m' ,m,'2 

- (mil M«(7p->A)| A,m') (m2| Mj^Hl m'^)] {m^m'^l 4/^)1^1^), 

where mI% and M^i are the elastic and charge-exchange {n^n — )■ tt p here) irNi amplitudes, 
respectively. The relative sign "-" between two terms in Eq. (fT7|) arises from isospin anti- 
symmetry of the DWF with respect to the nucleons. Further, we rewrite the denominator 
k'^-H^+iO in Eq. ([16]) as 



k''-lj'+iO = 2W2{E -E'+tO), E=^ki+m\ E' = k'^^ + tji^ (18) 

{k2 = 1^21, ^2 = 1^2!), where W2 is the effective mass of the rescattering 7rA^2 system, and 
E {E') is the total energy of the final (intermediate) nucleon in the 7[N2 rest frame. In a 
way similar to Subsection III CI we split the amplitude M^^^ into "on-shell" and "off-shell" 
parts, and obtain 

^(i)^^(i),on^^(i),o//^ /^^^ [A{^hu{p'2) + Iu)+B{^k2w{p'2) + Q], (19) 

fdk^uip^ fdk^wip^ 
" / nE' E'-E' f nE' E'-E' 



Here: dVt' = dz'dip' is the element of solid angle of relative motion in the intermediate vrA^ 
system; the factor A{B) is given by the r.h.s. of Eq. f lT7|) after the replacement \E'rf(p2) — )■ 
Su{Sw) (see Appendix IV Bt Eq. and is calculated with the on-shell intermediate pion 

and nucleon. The "off-shell" part M^'"^^ of the amplitude mP (pSD is given by the 
terms, containing the integrals lu and The vrA^-scattering amplitude is described in 
Appendix IV F[ 



III. RESULTS 



A. Comparison with the experiment 

We present herein the results of calculations and comparison with the experimental data 
on the differential cross sections da^did) / dVL, where VL and 6 are solid and polar angles of 
outgoing TT^'s in the laboratory frame, respectively, with z-axis along the photon beam. The 
results are given in Fig. |5]for a number of the photon energies E^. Calculations were done 



with DWF of the CD-Bonn potential (full model) 29| . The filled circles denote the data 
from the bubble chamber experiment at DESY [33|. 

The dotted curves show the results obtained with the lA amplitude [Fig. lH^a)]. It is 
known that the lA cross section (T{'yd-^TT^pp) can be expressed in the closure approxima- 
tion [34] through the cross section a{'~fn ^ n^p) and Pauli correction factor, which comes 
from the cross term of the amplitudes Ma^^ and M^^^ (|3]). It reads 



da , , , da 



dVt^' "^"^^ dVt 



K\ 



ild^7r-pp) = —i^n^7r-p) 1 - F^iA) + [■ ■ ■]=^^== F^iA) . (20) 



LP+ K 



Here: [■■■] is the Pauli factor, -^^(A) is the spherical form factor of the deuteron (we 
neglect the contribution of the quadrupole form factor), and A = is 3- momentum 



transfer; |Lp and \K\'^ are non-spin flip and spin flip 7n— J-vr p amplitudes, respectively (see 
Appendix IV Dj Eq. f l42|) ) squared and averaged over the photon polarization. For zero- angle 



{6= 0) pions, the non-spin flip term |Lp = 0. Then at A— t-O, we have ^^(A)— ^1 and the 
Pauh factor [■ ■ ■ ] ^ 2/3 in Eq. (EO]). At A ^ oo, we have ^^(A) ^ and [■ ■ ■ ] ^ 1. The 
momentum transfer A increases together with the laboratory angle 6. Thus, the spectra on 
Fig. [2] should be partly suppressed at small angles ^ ~ as compared with da/dQ{'yn-^TT~p). 
Fig. El^a) shows two different results for da^d{0)/d^l at E^ = 500 MeV: the dotted curve 



10 




FIG. 2: The differential cross section da/dD. of the reaction 'yd^TT~pp in the laboratory frame at 
different values of the photon laboratory energy <1900 MeV; 9 is the polar angle of the outgoing 
vr^. Dotted curves show the contributions from the lA amplitude Ma [Fig. H^a)]. Successive 
addition of the NN-FSI [Fig. ^h)] and vrA^-FSI [Fig. ^c)] amplitudes leads to dashed and sohd 
curves, respectively. The filled circles are the data from Ref. (33| . 
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FIG. 3: The differential cross section of the reaction jd^TT^pp in the laboratory frame at = 
500 MeV. (a): the dotted curve is the lA contribution, i.e., the same as in Fig. [21 the dashed one 
is the contribution from |Mi"'^^p + |M^^"'^^p, i.e., without the cross term, (b): the solid and dashed 
curves means the same as in Fig. [21 the dotted curve is obtained with the lA-term and ^-wave 



part of NN-FSI. The data are from Ref. 



33| 



represents the contribution from the lA amplitude squared |Mi^'' +M^^^^ p and the dashed 
one shows the contribution from |Mi^'' p+|M^^^'' p, i.e., without the cross term. The difference 
of the curves in Fig. [3]^a), i.e., the Pauli effect, at small angles is clearly seen. 

The dashed curves in Fig. |2] are the contribution of lA- and NN-FSI terms Ma + 
[Fig. [Tt^a),(b)]. The sohd curves show the results obtained with the full amplitude Ma+ Mb+ 
Mc [Fig. [I]^a),(b),(c)], including IA-, NN-, and vrN-FSI terms. Fig. [2] shows a sizeable FSI 
effect at small angles 6 < 30°, and it mainly comes from NN-FSI (the difference between 
dotted and dashed curves). Comparing dashed and solid curves, one finds that vrN-FSI 
affects the results very slightly. Note that at the energies = 300 — 500 MeV, the effective 
masses of the final np states predominantly lie in the A(1232)3/2'^ region. Thus, the plots 
at = 370 and 500 MeV of Fig. |2] show that the role of vrN-rescattering even in the 
A(1232)3/2+ region is very small. 
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Fig. H] demonstrates a reasonable description of the data [33| on da^did) I dVL. These data 



are also confirmed by recent results from the GDH experiments [35| in Mainz. Note that the 
data are absent at small angles 9 < 30°, where the FSI effects are sizeable. This is also the 
region of the most pronounced disagreements between the theoretical predictions of different 
authors [35]. 

The role of FSI is shown in more detail at = 500 MeV in Fig. [3]^b). Here, the dashed 
curve is the result obtained including the lA-term and the S-wave part of NN-FSI. The 
dotted and solid curves mean the same as in Fig. [2l Thus we see that at small angles, the 
S-wave part of NN-FSI dominates the FSI contribution. 

At large angles, the FSI effects are more significant as the photon energy increases. It 
is evident from the plots at > 1050 MeV in Fig. [21 Our interpretation is that at both 
high energies and at large angles, the role of configurations with fast final protons increases. 
For these configurations, the lA amplitude is suppressed by the deuteron wave function in 
comparison to the rescattering terms. These kinematical regions were considered in more 
detail in Ref. |l7|. 



B. Extraction of the jn^n p cross sections from the data 

The data on the deuteron target does not provide direct information on the differential 
cross section da/ dQ{'~fn ^ tt^ p) , because of the 'yd ^ n^pp squared amplitude term |M-ydP, 
where M^a = Ma + Mj, + Mc [Fig. [1] and can not be expressed directly through the term 
|My„p. Let us neglect for the time the FSI amplitudes, [Fig. [II^b),(c)] and let the final 
proton with momentum Pi(P2) t>e fast (slow) in the laboratory system and denoted by 
Pi{p2)- Then, the lA diagram Ma^^ with a slow proton pi emerging from the deuteron 
vertex dominates, Afi^^ is suppressed, and M^^ ~ Ma^\ This approximation corresponds to 
the "quasi-free" (QF) process on the neutron. In this case, one can relate the differential 
cross section da/dQi{'~fn — > n^p) on neutron with that on the deuteron target as follows. 
(Hereafter, Qi is the solid angle of relative motion in the n^pi pair.) From Eq. @, we get 

|Mi')|2 = 4m|M^J.^|2(27r)V(P2), {27rfp{p)=u^{p) + w^{p), jp{p)dp=l, (21) 

where p{p) is the momentum distribution in the deuteron. Making use of Eqs. ([2]) and 
fl2T|) . and multiplying by a factor of 2 (we include also the configuration when slow and fast 
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FIG. 4: The correction factor R, defined by Eq. (j27p . where 9i is the polar angle of the outgoing 
7r~ in the rest frame of the pair vr^+fast proton. The kinematical cut (j29p is applied. The solid 
(dashed) curves are obtained with both ttN- and NN-FSI (only NN-FSI), taken into account. 
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protons are replaced, and the amplitude Mi^-* dominates), we obtain 



dp^dQi dQi E^^'""'' ' ' E2 

(see, for example, Refs. jz-Ql). Here: Ei^ is the photon energy in the rest frame of the virtual 
neutron with momentum p' [Fig. [U^a)]; the factor E'^jE^ is the ratio of photon fluxes in 
^d and reactions; Q2 is the laboratory polar angle of final slow proton p2- Hereafter, 
we use the notation rfcr^^/rfpg'^^i) where index "i" specifies the ^d^-n'pp amplitude M^*^, 
namely M^^f = Ma^^ and M^^= Ma- The notation da^dl dp2dVti (without index) represents 
the differential cross section, calculated according to Eqs. (j2]) with full amplitude M^d = 
Ma + Mb + Mc. Let us rewrite Eqs. ([22]) in the form 

da^d , . do^n (lA) (full) 

- ^(P2) r r = rp rpsj, Vp = — — , rpgj = j-^, (23) 



dp^dn, ^"^^ dn^' 'P'^M, ^Qp^, 'i-M ^^p^y 

where for short, we use the notations (full) = da^dl dp2dVLi and {i) = dal^^/ dp2dVLi for 
i = QF and lA. Eqs. ( |23il enable one to extract the differential cross section da.yn/dQi on 
neutron from da-yd/dp2dQi, making use of the factors n{p2) and r. Here: the factor ^(pg), 
defined in Eqs. fl22|) . takes into account the distribution function p{p2) and Fermi-motion of 
neutron in the deuteron; r = Vprpgj is the correction coefficient, written as the product of 
two factors of different nature. The factor rp takes into account the difference of lA and 
QF approximations. Formally, we call it "Pauli correction" factor, since the lA amplitude 
= Ma^'^ + Ma^^ is antisymmetric over the final nucleons. However, the factors Vp in 
Eqs. ( |23l) and the expression in square brackets [■ ■ ■ ] in Eq. (l20l) are not identical. The 
factor Vpgj in Eqs. ( 123|) is the correction for "pure" FSI effect. 

Generally for a given photon energy E^, the cross section da^d/ dp2dVLi f l23|) with unpo- 
larized particles and the factor r depend on p2, 6*2, 9i, and y^i (4 variables), where 61 and 
ifi are the polar and azimuthal angles of relative motion in the final vr^pi pair. To simplify 
the analysis, we integrate the differential cross section on deuteron over P2 in a small region 
P2 < Pmax and average over ipi. Then, we define 

where the index "i" was introduced above (after Eqs. ( !22|) ). The cross section ( !24l) depends 



on Ery and 9i. We calculate the same integral from the r.h.s of Eqs. ( 1221) . Then, we take 
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the cross section a^n/dfli out of the integral / dp2, assuming n(p2) to be a sharper function. 
Thus, making use of Eqs. ( l22|) - flM|) . we obtain 



(25) 



where da^n/dVLi is averaged over the energy E'^ in some region E'^ ~ E^. The value c = 
c{pmax) can be called the "effective number" of neutrons with momenta p < Pmax in the 
deuteron. Under the restriction IP2I < Pmax in the integral for c ( 125|1 . we get 



Pmax 



cipmax) = J p{p)pdp^l at pmax^ 00. (26) 



A number of values of c{p) are given in the Table |T] for two versions of CD-Bonn DWF [2^. 
TABLE I: "Effective number" of neutrons with momenta p < Pmax in the deuteron. 



Pmax (MeV/c) 


50 


100 


200 


300 


Ref. 


c{Pmax) 


0.335 


0.719 


0.941 


0.981 


(full model) [29] 


c{prnax) 


0.326 


0.704 


0.932 


0.978 


(energy-independent) [29] 



Further, we rewrite Eqs. (]25ll in the form 

-^iE„e,)-cR^, R-RpR^si, ^P-Jqp-y ^^si - (27) 

Here: (i) = dal^^/dQi {i = QF and lA) and (full)= da^d/dQi (the definitions are different 
from those in Eqs. ( |23ll ): the factors R, Rp, and Rpsi similar to r, rp, and rp^j, 
respectively, but defined as the ratios of the "averaged" cross sections dal^^/dVti. 
Finally, we replace da^d/dVLi in Eqs. (1271) by the '-fd^ir^pp data and obtain 

-^{E,, e,) = c'\pmax) R~\E„ e,) -^{E,, e,), (28) 

where da^^^ /dVLi is the neutron cross section, extracted from the deuteron data da^^^ /dQi. 
Since the factor R =(full)/(QF) is the ratio of the calculated cross sections, we assume 
that (full) = da^Jli"' / dVLi = da^^f/dQi. The factor R in Eq. (EEl) is the function of the 
photon laboratory energy E^ and pion angle 6i in the Tc^pi frame, but also depends on the 
kinematical cuts applied. The value E^ in Eq. (!28|l is some "effective" value of the energy 
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E'^ = Ej{l + P cos 62) in the range E^{1±P). Limiting the momentum p2 to small values, we 
have /3 ^ 1 and E^ ^ E^. This approximation also improves, since p{p2) peaks at p2 = 0, 
where E'^ = E^. 

Eq. fl25]) is implied to be self-consistent, i.e., the '-fn-^ir^p amplitude, extracted from 
the da^^J' / dQi is the same as that used in calculations of the correction factor R. Then, 
the following iterations are proposed. The 1st step: one obtains the cross section da^^/dQi 
from Eq. fl28|) at R = 1 (no corrections), making use of the coefficient c{pmax), and extracts 
the 7n amplitude M^*^ (0th approximation). The next step: one calculates the factor R 
defined in Eqs. (^7^, making use of the amplitude M^*^ for the calculations of the cross 
sections dal^^/dQi. Then, we repeat the procedure of the previous step with new value of R, 
and obtain the amplitude M^]^ in the 1st approximation. The procedure can be continued. 
If the correction is small, i.e., R ^ 1 {\R — 1| ^ 1), then M^)^ is a good approximation for 
the corrected 'yn^n^p amplitude. Since, there are regions, where i? ~ 1 the FSI effects are 
insignificant and the preliminary analysis of the R factor is important for the procedure of 
the extraction of the •jn^n^p amplitudes. 



C. Numerical results for the R factor 



We present the results, obtained with the model discussed above, for the correction factor 
R, defined in Eqs. ( 1271) . The results depend on the kinematical cuts. We use cuts, similar 



to those applied to the CLAS data events 



and select configurations with 



IP2I < 200 MeV"/c < IpiI, (29) 

where Pi(P2) is the 3-momentum of fast (slow) final proton in the laboratory system. The 
results are given in Fig. [Has functions of the photon laboratory energy E^ and 61, where 
61 is the polar angle of outgoing vr" in the 7r~pi rest frame with z-axis directed along the 
photon momentum. 

The solid curves show the results for R, where the differential cross section (full) in 
Eqs. f l25|) takes into account the full amplitude Ma+ Mi,+ Mc [Fig. [Ha),(b),(c)]. The dashed 
curves were calculated, excluding the ttN-FSI contribution from the (full) cross section. The 
main features of the results in Fig. H] are 

1. a sizeable effect is observed in the region close to 6*1 = 0, which narrows as the energy 
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increases; 

2. the correction factor R is close to 1 (small effect) in the larger angular region. 

Since R consist of two factors Rp and Rpsi^ ^^^^ present them separately in Figs. El^a) 
and \5lc) for = 1000 and 2000 MeV, respectively. Here: dotted, dashed, and solid curves 
show the values of Rp, Rpsi, and R, respectively; the factor Rpsi was calculated with the 
full amplitude Ma+ Mf,+ Mc [Fig. [Tl^a),(b),(c)] taken into account. We find that Rp 7^ 1 at 
small angles, i.e., the factor Rp in addition to the pure FSI factor Rpsi also contributes to 
the total correction factor R. 

This can be naturally understood. Since Rp is the correction for the 2nd ( "suppressed" ) 
lA amplitude Ma^\ one should expect Ma^^ ~ Ma^^ and Rp 7^ 1 at ~ The probability 
of such configuration increases at ^i— t-O. It is clear that the possibility of the configuration 
Pi ~ P2 and the value of Rp should be rather sensitive to kinematical cuts. 

The dominant role of the S'-wave NN rescattering in the FSI effect was marked in Sub- 
section nil A[ This contribution to the factor R is presented in Figs. |5]^b) and [5]^d) for 
Ej = 1000 and 2000 MeV, respectively. Here, solid curves mean the same as in Fig. HI i.e. 
the total results; the dashed curves show the values R, where Rpsi takes into account only 
the correction from the S'-wave part of NN-FSI. Comparing the solid and dashed curves, we 
see that the FSI effect mostly comes from the S-wave part of pp-FSl. Note that the S-wave 
pp amplitude and the total elastic pp cross section aeiipp) sharply peak near the threshold 
at the relative momentum p^ ^ 23 MeV/c^. Thus, the S-wave NN-FSI effect should be 
important in some region ~ i-^-, at small angles as mentioned above and is evident 
from Figs. [5]^b) and [5](d). Obviously, the result is sensitive to the kinematical cuts. 



D. Factor R and Glauber approximation 

Now consider the region of large angles ^1, where FSI effects are small (i? ~ 1). In 
this case, we have the rescattering of fast pion and nucleon on the slow nucleon-spectator 
with small momentum transfer. Then, we may estimate the FSI amplitudes in the Glauber 



approach [36|, if the laboratory momentum of the rescattered particle ^ p (typical value 
in deuteron). For AW-FSI, this condition gives sin^^i ^ pWi/mE^, where Wi is the tt~pi 
effective mass. Taking p = 150 MeV/c, we get > 15.4°(10°) for E^ = 1000 (2000) MeV. 
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20 40 60 20 40 60 

(deg) 1^, (deg) 




20 40 60 20 40 60 



1^, (deg) 1^, (deg) 

FIG. 5: The correction factors at = 1000 MeV [(a),(b)] and = 2000 MeV [(c), (d)]. The 
sohd curves are the same as in Fig. [4l (a) and (c): the dashed (dotted) curves are the results for 
the factor RpsiiRp)i defined in Eq. (I27|) . (b) and (d): the dashed curves show the factor R, when 
Rpgj takes into account only the S-wave part of NN-FSI. 

As for the ttA^-FSI, we should also exclude some region close to 6i ~ 180°, where vr" is slow 
in the laboratory system. The high-energy A^A^-scattering amplitude can be written as 

M^^ = 2ip Wa'^j^ exp(6t) , (30) 
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where p, W, t, b, and are the relative momentum, AW effective mass, square of the 
4-momentum transfer, slope, and total AW cross section, respectively. The amplitude is 
assumed to be purely imaginary, and spin-flip term is neglected. Retaining only the S'-wave 
part of DWF, we obtain the lA- and NN-FSl amplitudes (M^ and M^) in the form 

= M(i) = (■■■) M, = {■■■)J, (31) 

/(Pp± 
(27r)2 

Here, the lA amplitude Ma^^ is equal to the 1st term in the r.h.s. of Eq. (jS]) with the 
replacement '^d{P2) ~^ u{p2)Su (see Eqs. (I35l)): the 2nd term {Ma ) of the lA amplitude 
is neglected; t = —b{p2±—p±)'^, where P2±ip±) is the transverse 2-momentum of slow final 
(intermediate) proton with Oz \\ p^ (fast-proton momentum). The factor exp(6t) is smooth 
in comparison with sharper DWF u{p±) in the integral J (!3T|) : thus, we neglect it for 
simplicity, i.e., calculate J ( !3T|) at 6 = 0. Considering the case of very slow proton-spectator 
with P2 ~ 0, we take u{p2) ~ u{0) for the lA term Ma in Eqs. ( 131]) . We also add the ttW-FSI 
amplitude Mc with the same assumptions as for the AW-FSI, i.e., Mc = —(1/4) (■ ■ • ) J . 



Finally, the FSI correction factor is R = \Ma+Mh+M^\^ /\Ma\'^, and with CD-Bonn DWF |29j, 
we obtain 

Here, we use some typical values cr^^ ~ 45 mb and cr^jy ~ 35 mb for the total cross sections 
at laboratory momentum piab ~ 1 — 1.5 GeV/c. For the integral J at 6 = in Eq. ( |32|) with 



CD-Bonn DWF |29j, one gets J = -(27r)"^^ a Inm^ in the notations of Eqs. ([36]). 

i 

Our Glauber-type calculations are extremely simplified in a number of ways and give only 
a qualitative estimation. Some predictions for the FSI corrections in the Glauber approach 
for TT^ photoproduction on light nuclei were done in Ref. js^]. The analysis 3^ of the 
reaction 'yd ix^pp at high energies of the photons, based on the approach of Ref. js^] , 
gave the Glauber FSI correction of the order of 20%. Similar values 15%-30% for this 



effect in the same approach were obtained in Refs. [2J, |25| , while our estimation (1321) gave 
smaller value ~ 5%. To comment on this difference in the results, let us point out the 
difference of the approaches used. Here, we use the diagrammatic technique. The analyses of 



Refs. 
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25I . I37I . |38| are based on the approach which considers a semi-classical propagation 



of final particles in the nuclear matter. The applicability of the latter approach to the 
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deuteron case is rather questionable. Notice that our approximate estimation in terms of 
Glauber FSI correction gives results similar to that obtained with our full dynamical model 
at large angles, i.e., the solid curves in Fig. HJ are in a reasonable agreement with the value 
of from Eq. 

Thus, we obtain the following behavior of the correction factor R, for the reaction 771—)- 
7r~p, calculated from the reaction 'yd^TT~pp at high-energy photon beam with slow proton- 
spectator. A sizeable effect 7^ 1 is observed in the relatively narrow region ~ 
dominated by the S'-wave part of NN-FSI with additional some contribution from the "Pauli 
effect" due to the "suppressed" lA diagram. Small but systematic effect |-R— 1| ^ 1 is found 
in the large angular region, where it can be estimated in the Glauber approach, except for 
narrow regions close to ~ or 61 ~ 180°. 

IV. CONCLUSION 

The incoherent pion photoproduction process 'yd — t- TT~pp was considered in a model 
containing the lA and FSI amplitudes. The NN- and vrA^-FSI were taken into account. 
The inputs to the model are the phenomenological — ttN, NN — )■ NN, and nN — )■ tcN 
amplitudes, the deuteron wave function, and the additional parameter for the off-shell 
behavior of the ^5*0 partial amplitude of pp-scattering. The Fermi- motion was also taken 
into account in the lA amplitudes as well as in the FSI {NN + vrA^) terms. 

The model reasonably describes the existing data on the differential cross section 
da/dQ{'~fd ^ Ti^pp). Sizeable FSI effects were observed at small laboratory angles 6 < 30° 
for outgoing pions, where the main part of the effect comes from the ^5*0 part of pp-FSl. 
In this angular range, the theoretical predictions of different authors reveal the most pro- 
nounced disagreements. Thus, future experiments on the reactions 'yd^nNN are welcome, 
especially at small angles 9 < 30°, where data are absent. 

The procedure to extract the differential cross section da/ dVL{'yn ^ ti~ p) on the neutron 
target from the deuteron data was derived in terms of the FSI correction factor r (1231) . To 
reduce the number of variables, we gave the results for the averaged correction factor R fl27j) . 
defined as the ratio of the differential cross sections da/ dQi{'~fd tt^ pp) , calculated with full 
amplitude as well as in the quasi-free-process approximation, where Qi is the solid angle of 
relative motion in the system 7r~+fast proton. Also the kinematical cuts with slow spectator 
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proton were used. The results show a sizeable FSI effect R ^ 1, predominantly coming from 
the ^5*0 part of pp-FSl, at the angular region close to ~ 0, and the region narrows with 
the increasing photon energy. In the wide angular range, the effect is small 1| ^ 1) 
and in agreement with the Glauber estimations. 

The more refined analysis requires the use of the factor r (125]) instead of the averaged 
one (R). Then, we deal with the ratio of multi-dimensional differential cross sections 
dal^^/dp2dQi, used in Eqs. fl23|) . Further, one should integrate dal^^/dp2dQi over the az- 
imuthal angle (pi in the ir^pi pair, since the differential cross section on the neutron in the 
unpolarized case has no azimuthal dependence; thus, the cross sections dal^^/ dp2dVLi turns 
out to be a function of 3 variables, i.e., p2, 6*2, and 9i (or cos6'2 and cos^^i). Thus, applying 
Eqs. f l2^ to extract the differential cross section da^n/dVli on the neutron, one needs data 
on the deuteron cross section da^d/dp2dQi binned in the variables p2, O2, and di, i.e., in the 
3-dimensional form. We plan to discuss this question in detail in the next publication. 
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V. APPENDIX 

A. Invariant amplitudes and phase space 

We use standard definitions and the cross section of the process a + b~^l +■ ■ -+71 reads 

CTn = InJ-' [ |M|2rfr„, dTr, = {27,f5^'\P, - Pf) TT TX^TT- (33) 

Here, M is the invariant amplitude; (ir„ is the element of the final n-particle phase space; 
Pi{Pf) is the total initial (final) 4- momentum; E^^ and Pj are the total energy and 3- 
momentum of the i-th final particle; J = AEarrib = iqaby/s is the flux factor, where Ea 
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{mb) is the total laboratory energy (mass) of the particle a{b), qab is the initial relative 
momentum and ^/s is the total CM energy; /„ = l/rii! ■ ■ -rifc! is the identity factor, where 
Hi is the number of particles of the i-th type (ni + - ■ ■ + nk = n). 



B. Deuteron vertex and wave function 

The deuteron vertex F^, used in Eq. (j3]), can be written in the form 

Up) = /^(ep) + - 

3m2 (34) 

gi = — —y/m{p^ + a'^)w{p), g2 = \/m{p'^+a'^)[\/2u{p)+w{p)]. 

Here, e is the deuteron polarization 4-vector; p = \p\ the relative 3-momentum of the 
nucleons; u{p) and w{p) are S- and D-wave parts of the deuteron wave function, respectively; 
= med, where is the deuteron binding energy. The DWF in p-representation reads 

{mi,Ti,m2,T2\ ^d{p)\md) = (ff^dip) V^E, 

(o-e) 1 (35) 

"ifdip) = u{p)Su + w{p)Sw, Su = —j=-, = -[(ere) - 3{ne){an)]. 

Here, n = p/p; e is the deuteron polarization 3-vector for a given spin state m^; and Xj 
are spin and isospin states of the i-th nucleon, and (p^ is its spinor and isospinor; y?^ = r20"2V?*, 
where a2 and r2 are spin and isospin Pauli matrices. We use the normalization 

I [dp Yl Ti, 1712, T2\ ^ d{p)\ md)\^ = J dp [u\p)+w\p)] = {271 f. 

mi,Ti,m2,T2 

For the DWF of the CD-Bonn potential, the functions u{p) and w{p) were parameterized 29 1 
in the form 

The parameters q, di, and mj are given in the Tables 11 (full model) and 13 (energy- 



independent model) of Ref. [29 1 



C. Invariant ^N^ttN amplitudes 

The general expression for the 'jN^ttN amplitude M^tv can be written as 



4 



M^N = u{p2)M'yNu{pi), M-yN = i ^ Alb^i, (37) 



1=1 
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where u{pi2) are the nucleon Dirac spinors {uu = 2m), Ai are the invariant amphtudes, Fj 
are the 4x4 matrices. Fj's can be taken in the form 

ri = #, T2 = {ep){qk) - {pq){ek), ^^^^ 

^3 = i{ek)-^{qk), T 4 = ^{ep) - ^{pq) , p = p^ + p^. 

Here, e is the photon polarization 4-vector; q, k, and pi^2 are 4-momenta of the photon, pion, 
and nucleons, respectively. One can write the amplitude M^n (!37I1 in CM frame as 

M^N = SirW (ptF(pi, F = ieFi + n2[(T{niX e)]F2 + ihi{n2e)Fs + in2{n2e)F4. (39) 

Here, e is the photon polarization 3- vector; q*{k*) are the photon (pion) CM 3-momenta; 

W is the total CM energy; Fi = Fi{W, z) are the CGLN [32] amplitudes, z= cos 9; ip^ are 

the Pauli spinors; rii = q*/q*, n2 = k*/k*, q* = \q*\, k* = \k*\; "hat" means the product 

k* 

with cr, i.e., e = (ere), etc. For unpolarized nucleons da/dQ{'~fN ttN) = Tr{FF~^}. 

Equating Eqs. (l37|l with Eqs. (139|) . one finds the relations between Aj's and Fj's, i.e., 

, _ F, + F2 p^S^^^p - _ 8nWN,N2 

' 2W ' ' N,N2 ' \q\\k\ 

^'-^mT' \q\\k\W^N2^'' \k\m.Nj'' 

_ F2+iy+Ai-(A:g)A34 , _ W^F^^W.F, 
A3-A4 + /I34, A- , A34- ^ , 



where W± = W ±m, Ni^2 = \fF^L'2?rrn^ and Ex^2 are total CM energies of the nucleons. 

The isospin structure of the amplitudes Aj(7A^— j-vTaA^) and contributions to the different 
charge channels read 

A^ = A'^^S^s + A^-^Ik r,] + AfK^ , (41) 
Ai^p^TT^p) = A^^ + Af\ Ai^p^TT+n) = V2{Af^ + 4"^), 



Ai {jn^7r^n) = A^^^ - ^ , Ai (772 ^n-p) = V2{AXI - A 



(0) .(-)^ 



The amplitudes Ai{'yN -^ttN) can be obtained from the CGLN 
ttN) through Eqs. (HOl) . We use the GW pion photoproduction amplitudes Fi 



32| amplitudes Fj(7A^— )• 



26| 



D. Matrix elements for jN^ttN 

The matrix element (m2| M-yn\ A, mi) in an arbitrary frame can be written in the form 



{m2\M^n\\mi) = NxN2{m2\L + t{K(T)\mi) {Ni = y/Ei+m). (42) 
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Making use of Eqs. ( I37|) - (l38|) . we obtain 

L = Ai{e[q x (xi-xa)]) - (er)(qr[xiX xa]) + [Aiqo - (gr)](e[xiX xa]), 
K= [AiCi + (gr)c3]e + {eS)q + (e52)xi + (e5i)x2 + A2C2(x2- xi), 

Xl,2 =Pl,2/(^l,2+m), ^^^^ 

S = Ai(xi+X2) + cgr, 5i,2 = [(gr) - Aigo]xi,2 + [(qxi_2) - go]^, 
Ci = go(l + X1X2), C2 = 2[{qpi){ek) - {qk){ep^)], C3 = 1 - (X1X2), 
r = A3k + Ai{pi + p2), r = A3k + A^ip^+p^). 

Here, Ai are the amplitudes in Eqs. fl57l) : e = e*^'^^ is the photon 3- vector, specified by spin 
state A; 5,^,^1,2(9,^,^1,2) the 4(3)-momenta, defined in Appendix IV CI We fix two 
possible photon states (A = 1, 2) by definition e\^^= 6ix, where e^^'^^ is the i-th component of 
e(^) {Oz II q). Thus, (e(i)e(2)) = and (e^q) = 0. 



E. Invariant A^A^ — > A^A^ amplitudes 

The NN-scattering matrix depends on 5 independent spin amplitudes, and different 
choices can be found in Refs. 0, |4^. In the NN rest frame, the N[N2 — ?■ N1N2 matrix 
element can be written in the form (m^,m2| Mnn\ 'mf,m2') = 87iW{Fi\in), where W is the 
AW effective mass, and 
4 

{Fnn) = fiiftQif'i)iftQif2) + Miftnfi){ftv>2) + iftfi)iftnf'2)], (44) 

i=l 

where ip[ 2 (v^i 2) the Pauli spinors of the initial (final) nucleons, specified by spin states 



m 



^2 ('^12)- Here, we use the formalism of Ref. [40 



where /i, ■ ■ ■ /s are the independent 



spin amplitudes; Qi, - ■ - Qa are the 2x2 matrices, and 

Qi=/, Q2=n, Q3=m, Q4=/, n J^y^*} m=ff^ i=f±J^. (45) 

p*' = Pi' — P2 {p* = Pi~ P2) is the initial(final) relative momentum. 

To apply Eq. fj44|) for calculation of the matrix elements (mi, m2| Mat at | itl[, m'2) in Eq. ([6]), 
one should transform the NN amplitude from the deuteron rest frame to the NN rest frame. 
The possible way is to transform the nucleon Dirac spinors to the NN rest frame, and find 
the corresponding unitary transformation of spinors in Eq. fl44p . i.e., 

U = N-\L + iK(t), N=^/\VfT\K^, (46) 
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where is any of V^'i 2- The result is 

L = ao + &x, X = a + feox + [6 x x], x 



E' + m' 

ao = C1C2, a = (— S1S2, C1S2, S1C2), &o = ~2;jy^SiC2, & = — a;jY7v(ciS2, S1S2, C1C2), 



(47) 



9, 



■ '■PnN ^NN ■ '^NN ■ "NN 

Si = sm ^ , ci = cos , S2 = sm— — , C2 = sin— z 



2 ' - 2 ' 2 ' - 2 ' ^iviv+W^" 

Here: and p are the total energy and 3-momentum of a given nucleon in the deuteron 
rest frame, i.e., p = Pi,2,p'i,2 Pig- H^b)]; -^aw, Patat, and v^atat are the total energy, 
3-momentum, polar and azimuthal angles of the outgoing AW system in the deuteron rest 
frame, respectively. Finally, for the AW matrix elements in Eq. (jG]), we obtain 



4 

(Fnn) = Y,fi (mil f/+4f>{| m[){m2\ U^QiU^lm',) 
1=1 



(48) 



One can rewrite the products UQU' in the form UQU' = Vq + iVcr, making use of 
Eqs. ( H5|l - (H7|) . and calculate the factors (mj| ■ ■ ■ | m^) in Eqs. fHHj) (we ommit the details). 



The Hoshizaki 



40| amplitudes /i, ■ ■ - /s can be expressed through the helicity amphtudes 



Hi, - ■ ■ Hr, (the relations of iJ's to other representations 39 



40| can be found, for example. 



in Ref. [411), the results of GW NN partial- wave analysis [27 1. 



F. Invariant vrA^ — vrA^ amplitudes 

Calculating the ttN ^ttN matrix elements in arbitrary frame, we start from the invariant 
amplitude and write 

M^^ = U2{A + B^)ui = ^t^Vi- (49) 

Here, tti,2(v'i,2) are Dirac (Pauli) spinors; A and B are the invariant amplitudes; p = {po,p) = 
Pi, 2+ ^1,2 is the total 4-momentum; pi,2= (-£'1,2, Pi 2) (^1,2) are the 4-momenta of the initial 
and final nucleons (pions); □ is 2 x 2 matrix. Making use of Eq. f H9l) . we obtain 

D = N{L + iKa), N= ^{Ei+ m){E2+ m), 
L = A + Bpo- B{p (X1+X2)) + {Bpo- A)(xiX2), (50) 
K = B[px (X2-X1)] + {Bpo- A) [xi XX2], xi,2 = Pi,2/(^i,2+ m). 
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The matrix elements can be obtained from Eqs. flSUp . i.e., {m2\M^^\ mi) = (m2|[I]| mi). 

In the ttN rest frame, □ = 8ttW[F+ iG{[nixn2](T)], where F{G) is the standard non-flip 
(spin-flip) amplitude, W is the effective vrA^ mass, ni^2 = Pi 2/IPi 2I) Pi 2 nucleon CM 

3- momenta. Applying Eq. (jUj), one can relate the amplitudes A and B to F and G, i.e., 

A.,.W{^,1), B^,.{^-1), E^^E±rn, (51) 

where E is the nucleon total CM energy, z is the cosine of CM scattering an gle. We use the 
amplitudes F and G, based on the results of GW nN partial- wave analysis [28j] . 
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